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f— ^ ' Abstract. Deep two-color imaging of 40 clusters of galaxies with the NOT 

^S| , and the University of Hawaii 2.24-ni telescope is used to measure the weak 

gravitational shear acting on faint background galaxies. From this, maps 

5J^. of the projected cluster mass distribution are constructed, and the cluster 

^Q ' masses are measured within circular apertures of up to ~ 3/i~^ Mpc. The 

ly-v , results are used to derive the cluster mass function at z ~ 0.2. The average 

^SJ ' mass-to-light ratio of the clusters indicate a low-density Universe with fio — 

0.25. Additional mass concentrations are found along the lines of sight to the 
clusters, some of which are most likely associated with clusters at z > 0.5. 
In addition, the data are used to measure weak gravitational lensing due to 
large-scale structures along the line of sight. The future prospects for this 
^r^ , kind of research with the NOT are briefly discussed. 
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The weak gravitational lensing phenomenon has for a decade been used to probe the 
Q^. dark matter distribution in the Universe on scales ranging from galaxies to large-scale 

structures. Rapid progress has been made both observationally and theoretically, and 
many important developments have taken place only during the last 2-3 years. Much 
of the progress on the observational side is due to the development of large-format 
53 . (usually mosaic) CCD cameras which can cover fields of order 0.5 degrees on a side 

on 4-meter class telescopes. Such cameras are currently among the most heavily 
used instruments at leading observatories, and significant resources are being spent 
^\ ' on building even larger CCD mosaics. Weak lensing studies are frequently mentioned 

as an important part of the scientific motivation for these technical developments. 

On the theoretical side, many different inversion algorithms for reconstructing 
the mass distribution from measurements of weak gravitational shear have been de- 
veloped, e.g. to take into account strong lensing in the inner regions of the clusters. 
With the increase in the size of CCD cameras, weaker and weaker shear can in princi- 
ple be measured, provided that systematic effects can be kept under control. Recently 
developed algorithms can account for the effects of realistic PSFs and optimally weigh 
background galaxies of different fiuxes, shapes and sizes when estimating the weak 
shear. 

Why all the interest in weak lensing? For the first time, we are getting an 
unbiased (although still fairly noisy) look at the mass distribution in the Universe. 
We are now at a stage where measurements are getting sufficiently sensitive to map 
the projected mass distribution along any line of sight. In many ways, this is similar to 
what happened in the past when new parts of the electromagnetic spectrum became 
available for astrophysical studies, but the mass probed by lensing is a much more 
fundamental quantity than electromagnetic radiation emission and this makes weak 
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lensing measurements a very powerful tool for distinguishing between different models 
for structure formation. 



2. The cluster sample 

Our observational efforts at the NOT have been focused on measuring weak gravita- 
tional shear in cluster fields. Although many weak lensing studies over the past decade 
have targeted clusters of galaxies - see e.g. the recent reviews by Mellier (1999) and 
Bartelmann & Schneider (2000) and references therein - they have mostly been lim- 
ited to studies of a few clusters which are thought to be extremely massive, thereby 
guaranteeing a significant weak lensing signal. Until now, there has been no weak 
lensing study of a large, well-defined sample of objects. The sample described here 
was selected from the lists of X-ray luminous clusters published by Briel & Henry 
(1993), Ebeling et al. (1996; XBACS) and Ebcling et al. (1998; ECS). The clusters 
are in the redshift range 0.15 < z < 0.35 and almost all have X-ray luminosities 
Lx,o.i-2.4koV > 10''''' erg s^^ (for h ~ 0.5, go ~ 0.5). Most of the clusters are targets 
for recent or planned X-ray observations with Chandra or XMM-Newton. 



3. Observations and reduction 

The NOT observations were made with ALFOSC in imaging mode, and 1.5h of expo- 
sure was typically obtained for each cluster in each of the V- and I-bands. Additional 
observations were made with the University of Hawaii 2.24-m telescope at Mauna 
Kca Observatory with either a 2k Tektronix CCD camera or the UH8K CCD mosaic 
(Luppino et al. 1995). The UH8K pixels were 2x2 re-binned to effectively make a 4k 
detector with a 19' field. To avoid introducing artificial shape distortions mimicking 
the effect of weak lensing, the individual exposures had to be accurately aligned to 
produce a deep combined exposure of each field. For the 2k data, catalogs of stellar 
positions was generated for all exposures, and an exposure with minimal extinction 
was chosen to be the astrometric and photometric reference for the other exposures. 
A least-squares solution was then calculated for a lower-order polynomial astromet- 
ric transformation of each exposure into the reference, and the images were median 
combined into a final image. To eliminate artificial galaxy shape distortions due to 
PSF anisotropics, a spatially variable PSF model was calculated, and the image was 
convolved with the model, rotated by 7r/2, to yield a smoothed image with a nearly 
circular PSF. The PSF model, along with the combined image and the smoothed 
image, was used to estimate the gravitational shear acting on a sample of faint back- 
ground galaxies in the field with optimized weighting of the galaxies as a function of 
their fiuxes, sizes and shapes (for details on the shear estimator, see Kaiser 2000). 

The 8k data were reduced using a proceduree very similar to the one described by 
Kaiser et al. (1999). An astrometric solution was found which accounts for the CCD 
chip layout in the camera and field distortions (which were quite small), using the 
USNO-A2.0 catalog as a reference. The PSF is known to vary discontinuously across 
chip boundaries in the UH8K camera, so the individual chip frames were convolved 
with a re-circularizing filter generated as described above before combining the frames. 
A set of 'raw' and 're-circularized' combined images were then made, and the shear 
was estimated in essentially the same way as for the 2k data. 
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KS93 ^(r) RECONSTRUCTION OF rxj1532 
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Figure 1. The upper left panel shows an I-band image of the cluster 
RXJ1532.9+3021 (z^O.345), based on 3 hours of exposure with ALFOSC in May 
1999. The upper right panel shows the projected mass density in the field, calculated 
from weak lensing measurements. The lower left panel shows the predicted mass den- 
sity, based on photometry of early-type cluster galaxies, and the lower right panel 
shows the number density of galaxies in the field. The contour plots have all been 
smoothed with a Gaussian of scale 36", and the average value has been subtracted 
from the two lower plots. 
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4. Weak Lensing Analysis 

4.1. Predicted and reconstructed surface density 

Using local measurements of the spectral energy distributions (SEDs) of galaxies of 
different types (Coleman et al. 1980), and assuming negligible galaxy color evolution, 
it is possible to predict the broadband colors of galaxies as a function of redshift. 
Lubin (1996) shows that the colors of early-type galaxies at redshifts z < 0.4 are 
well represented by this no-evolution prediction. In this redshift range, early-type 
galaxies can be separated by their V-I colors from spirals and irregulars, and early- 
type galaxies in rich clusters form a tight, well-defined sequence in color-magnitude 
space. It is thus fairly easy to identify and isolate the early-type cluster galaxies and 
other early-type galaxies at moderate redshifts. Assuming that the mass traces the 
early-type galaxies with a constant mass-to- (early type galaxy) light ratio, a map 
of the predicted mass distribution in the field can be generated and compared to 
the mass distribution derived from weak lensing. The projected surface density was 
reconstructed from the shear measurements using the algorithm of Kaiser & Squires 
(1993), which is fast to implement and has well-understood noise properties. An 
example of the results for the cluster RX J 1532. 9-1-3021 is shown in figure |[ 

4.2. Light-mass cross-correlation 

To investigate the relative distribution of mass and galaxy light in the field, the two- 
dimensional cross-correlation between the measured and predicted mass distribution 
was calculated. This was compared to the auto-correlation of the predicted mass 
distribution from galaxy photometry. The results of this procedure when applied to 
wide-field data is particularly interesting, since they probe a wide range of physical 
radii in the clusters. The clusters display a range of properties: In some clusters, the 
dark matter appears to be significantly more concentrated than the early-type galaxy 
light, in others the two components trace each other within the errors, and in some 
clusters the mass is significantly more extended than the light. The differences appear 
to be related to the dynamical state of the clusters. The clusters which appear to 
be completely relaxed systems have the more concentrated dark matter distributions, 
whereas the clusters with the more extended dark halos appear to still be in the later 
stages of their formation process. 

The rest-frame B-band mass-to- light ratios in the cluster fields were calculated by 
taking the discrete Fourier transform of the predicted and measured surface density 
images and comparing the amplitudes of the Fourier modes. In figure 0, the resulting 
values for M/L-q are plotted against the physical scale corresponding to the longest 
wavelength measured. Also plotted are a number of other M/Lb measurements taken 
from the literature (Bahcall et al. 1995). The average value for the observed clusters 
is M/Lb = 377 ± 17 in solar units, which coupled with local measurements of the 
B-band luminosity density (Loveday et al. 1992) indicates a value for the density 
parameter of fio — 0.25. 

4.3. Comparison with other mass measurements 

Many of the clusters in the sample studied here have X-ray temperature measure- 
ments and/or galaxy velocity dispersion measurements available. By fitting a singular 
isothermal sphere model to the measured tangential shear profiles around the clusters, 
the weak lensing measurements can be compared to velocity dispersion measurements, 
as shown in figure |3|. The agreement is generally fairly good, with a tendency for the 
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Figure 2. (Rest-frame) B-band mass-to-light ratio vs. physical scale. The M/Lb 
expected for fl — 1 and fl — 0.3 are shown as lines, along with empirical Rts to M/Lb 
(assumed to increase proportional to R up to a scale oflOOh^^ kpc) for spirals and 
ellipticals. The plotted points for rich clusters, Morgan groups, Hickson groups, CfA 
groups, spirals and ellipticals are median values of these samples. The circles are 
values for the clusters studied in this work, using the mass to total light ratio. The 
figure is in part based on a similar plot by Bahcall et al. (1995). See this paper for 
further references. 
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Figure 3. Spectroscopically measured velocity dispersion (t„ vs. the dark matter 
velocity dispersion ctdm , inferred from weak lensing measurements. 

spectroscopically measured velocity dispersion to be marginally larger than the dark 
matter velocity dispersion (derived from the weak lensing measurements) for about 
half of the clusters. This is not unexpected, since the spectroscopically measured 
velocity dispersions may be biased upwards due to line-of-sight contamination by 
galaxies in the infall region outside the virial radius of the clusters. 

4.4. The mass function of galaxy clusters 

From the list of observed clusters, a volume-limited BCS subsample with the selection 
criteria 0.1 < z < 0.25 and Lx,o.i-2.4kcV > 10^^ erg s^^ was selected. Using the singu- 
lar isothermal sphere model fit, the masses contained within the inner 0.5h~^ Mpc of 
the clusters were calculated, and a cumulative cluster mass function was constructed, 
as shown in figure IJ. This is the first mass function of galaxy clusters which is based 
on gravitational lensing measurements. Further work is currently in progress to com- 
pare the results to predictions from Press-Schechter theory and from recent large 
cosmological N-body simulations (e.g. Jenkins et al. 2000). 

4.5. Detection of additional mass concentrations along the line of sight 

In addition to the targeted clusters, some additional mass concentrations are seen 
in the reconstructed density maps of some cluster fields. These mass concentrations 
tend to be associated with concentrations of red {V — I > 2.2), faint galaxies, which 
indicates that they are distant clusters (z > 0.5). One distant cluster candidate, 
situated in the field of Abell 1705, even shows a strongly lensed blue arc surrounding 
what appears to be the central cluster galaxy (see figure pi). Thus, the weak lensing 
measurements can be used to discover new clusters in the observed fields. The avail- 
ability of color information is a great help in this respect: What initially appeared 
as a 'dark cluster' in the field of Abell 1722, similar to the one reported by Erben et 
al. (2000) , turned out to be associated with a concentration of red galaxies which did 
not stand out from the background/foreground in data from a single passband only. 
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Figure 4. Cluster mass function. The quantity A/(0.5) is the mass within a cluster 
radius ofO.Bh^^ Mpc in three-dimensional space. 




Figure 5. A blue arc curving around the central galaxy of a previously unknown, 
distant cluster. 
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4.6. Cosmic shear in cluster fields 

Several groups have recently reported the detection of weak cosmic shear from large- 
scale structures along the line of sight (Bacon et al. 2000; Kaiser, Wilson & Luppino 
2000; van Waerbeke et al. 2000, Wittman et al. 2000). These studies are generally 
based on observations of 'blank fields' not containing any previously known clusters. 
The presence of the massive clusters in the imaging data described here may bias 
the cosmic shear measurement somewhat, but the shear on different sides of the 
cluster (which is generally situated near the center of each field) will mostly cancel 
out, and measurements of the net shear across the image will thus only be weakly 
affected by the presence of the cluster. Using 12 cluster fields observed in the I-band 
with the UH8K camera, the cosmic shear variance on scales of 20' was measured to 
be (7^)// = 1.24 ± 0.49 x 10"'*, consistent with the results of Wittman et al., but 
marginally larger than the Kaiser et al. result. Using a cross-correlation of shear 
measurements from I-band data and shallower V-band data in 7 UH8K fields (which 
will cancel out systematic effects which are not common to both passbands) , an upper 
limit of (7^)/v < 0.50 x 10^** was derived, more in line with the values found by Kaiser 
et al. 



5. Future prospects: What NOT can (and can't) do 

Although the Mosaic/FRED combination would still be extremely competitive if both 
instruments were already available, the NOT will within a couple of years become less 
attractive for weak lensing studies requiring very wide fields (e.g. measurements of 
cosmic shear and weak galaxy-galaxy lensing) , since the field sizes of competing facil- 
ities (such as Subaru Suprime-Cam, CFHT MegaPrime, MMT Megacam, VST and 
VISTA) are rapidly increasing. In the more distant future, wide-field weak lens- 
ing observations will most likely be c arried out using dedicated facilit ies such as the 



proposed Dark Matter Telescope (see http: //www.dmtelescope . org) or an array of 



smaller telescopes (Kaiser, Tonry & Luppino 2000). Until then, the NOT (with Mo- 
saic/FRED) will likely remain competitive for weak lensing observations that target a 
fairly large number of objects along widely separated lines of sight, such as studies of 
known galaxy clusters and superclusters. The main advantages of the NOT will still 
be the excellent intrinsic seeing at the telescope site (work to further improve image 
quality should have very high priority !!) and the possibility for relatively extensive 
telescope time allocations. 
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